function r = quants(r0, v0)
global g_delta_r;
global g_delta_phi;
global g_max_radius;
global g_survey_num;
global g_period;

delta_r_num = floor(g_max_radius / g_delta_r);
delta_phi_num = floor(2 * pi / g_delta_phi);

% Шум
%r = zeros(g_survey_num, delta_r_num * delta_phi_num);
r = floor(random('Uniform', 0, 1.05, g_survey_num, ... 
    delta_r_num * delta_phi_num));

% Истинные отметки
for survey = 1:g_survey_num
    ri = r0 + (survey - 1) * v0 * g_period;
    index_r = floor(abs(ri) / g_delta_r);
    index_phi = ceil(angle(ri) / g_delta_phi) + 1;
    index_total = index_r * delta_phi_num + index_phi;
    if (index_total <= delta_r_num * delta_phi_num && index_total > 0)
        r(survey, index_total) = 1;
    end
end

end